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Abstract 

QCD corrections to order a s in the soft-gluon approximation to an- 
gular distributions of decay charged leptons in the process e + e~ — ► it, 
followed by semileptonic decay of t or t, are obtained in the e + e~ 
centre-of-mass frame. As compared to distributions in the top rest 
frame, these have the advantage that they would allow direct com- 
parison with experiment without the need to reconstruct the top rest 
frame. The results also do not depend on the choice of a spin quan- 
tization axis for t or t. Analytic expression for the triple distribution 
in the polar angle of t and polar and azimuthal angles of the lepton 
is obtained. Analytic expression is also derived for the distribution 
in the charged-lepton polar angle. Numerical values are discussed for 
^/s = 400 GeV, 800 GeV and 1500 GeV. 

PACS: 12.38 -t, 14.65.Ha, 13.65. +i 

Keywords: Top quark, electron-positron annihilation, QCD correc- 
tions 



1 Introduction 

The discovery [0] of a heavy top quark, with a mass of about 174 GeV which 
is close to the electroweak symmetry breaking scale, raises the interesting 
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possibility that the study of its properties will provide hints to the mecha- 
nism of symmetry breaking. While most of the gross properties of the top 
quark will be investigated at the Tevatron and LHC, more accurate determi- 
nation of its couplings will have to await the construction of a linear e + e~ 
collider. The prospects of the construction of such a linear collider, which 
will provide detailed information also on the W ± , Z and Higgs, are under 
intense discussion currently, and it is very important at the present time to 
focus on the details of the physics issues ( |2j and references therein). 

In this context, top polarization is of great interest. There has been a lot 
of work on production of polarized top quarks in the standard model (SM) in 
hadron || collisions, and e + e~ collisions in the continuum [4J], as well as at 
the threshold []5J. A comparison of the theoretical predictions for single-top 
polarization as well as ti spin correlations with experiment can provide a 
verification of SM couplings and QCD corrections, or give clues to possible 
new physics beyond SM in the couplings of the top quark [6-18] (see |L9| for 
a review of CP violation in top physics). 

Undoubtedly, the study of the top polarization is possible because of its 
large mass, which ensures that the top decays fast enough for spin informa- 



tion not to be lost due to hadronization [20]. Thus, kinematic distributions 



of top decay products can be analysed to obtain polarization information. 
However, most studies on top polarization, particularly on QCD corrections, 
have presumed an accurate reconstruction of the spin, and are generally not 
concerned with decay distributions. Some works do consider decay distri- 



butions, but in the top rest frame JT8|, g^, 23|. In these cases, the onus 



of accurate determination of the spin quantization axis and that of the top 
energy-momentum and hence of the top rest frame, is largely left on the ex- 
perimentalists. The remaining papers which discuss decay distributions in 
the laboratory frame are restricted to energy distributions of decay products 
[[HJ, and do not calculate angular distributions. 

The choice of spin quantization axis has also been an issue of discussion 



[ 21] , |22|, and it has been remarked that a certain "off-diagonal" spin basis 
has certain advantages. 

In an alternative approach, adopted usually in the context of the study 
of top polarization arising from new physics, predictions are made directly 
for decay-lepton 0, §, g [H| 0, [12], [14], [15], |I§ or bottom-quark 0, 
distributions in the laboratory frame. Such an approach makes the issue of 
the choice of spin basis for the top quark superfluous. Moreover, if the study 
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is restricted to energy and polar angle distributions of top decay products, it 
even obviates the need for accurate determination of the energy or momentum 
direction of the top quark |TIJ. Even in case azimuthal distributions are 
studied, where an additional direction is needed to define the azimuthal angle, 
the knowledge of the magnitude of the top momentum is not required. It is 
sufficient if the direction is determined with reasonable accuracy. 

In this paper we shall be concerned with the laboratory-frame angular 
distribution of secondary leptons arising from the decay of the top quarks in 
e + e~ — > it in the context of QCD corrections to order a s . QCD corrections to 
top polarization in e + e~ — > it have been calculated earlier by many groups 
[23fl-|[30|. In cases when the corrections have been applied to the process 
including top decay, the discussion of angular distributions is restricted to 
the top rest frame f23fl . QCD corrections to the lepton energy distributions 



have been treated in the top rest frame in [T8j and in the lab. frame in 



[ p~3|j . This paper provides, for the first time, angular distribution in the e + e~ 
centre-of-mass (cm.) frame. As a first approach, this work is restricted, for 
simplicity, to the soft-gluon approximation (SGA). SGA has been found to 
give a satisfactory description of top polarization in single-top production 
and it is hoped that it will suffice to give a reasonable quantitative 
description. 

The study of the laboratory (lab.) frame angular distribution of sec- 
ondary leptons, besides admitting direct experimental observation, has an- 



other advantage. It has been found |L5], that the angular distribution is 



not altered, to first-order approximation, by modifications of the tbW decay 
vertex, provided the 6-quark mass is neglected. Thus, our result would hold 
to a high degree of accuracy even when 0(a s ) soft-gluon QCD corrections 
to top decay are included, since these can be represented by the same form 
factors [0 considered in |15|, |T6] . We do not, therefore, need to calculate 



these explicitly. It is sufficient to include 0(a s ) corrections to the ^yit and 
Zit vertices. This, of course, assumes that QCD corrections of the nonfac- 
torizable type f33j, where a virtual gluon is exchanged gluon between t (i) 
and b (b) from i (t) decay, can be neglected. We have assumed that these are 
negligible. 

The procedure adopted here is as follows. We make use of effective ■yii and 
Zit vertices derived in earlier works in the soft-gluon approximation, using 
an appropriate cut-off on the soft-gluon energy. In principle, these effec- 
tive vertices are obtained by suitably cancelling the infra-red divergences in 
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the virtual- gluon contribution to the differential cross section for e + e~ — > it 
against the real soft-gluon contribution to the differential cross section for 
e + e~ — > tig. For practical purposes, restricting to SGA, it is sufficient to 
modify the tree-level jit and Zti vertices suitably to produce the desired re- 
sult. Thus, assuming 0(a s ) effective SGA vertices, we have obtained helicity 
amplitudes for e + e~ — > it, and hence spin-density matrices for production. 
This implies an assumption that these effective vertices provide, in SGA, a 
correct approximate description of the off-diagonal density matrix elements 
as well as the diagonal ones entering the differential cross sections. Justifi- 
cation for this would need explicit calculation of hard-gluon effects, and is 
beyond the scope of this work. 

We have considered three possibilities, corresponding to the electron beam 
being unpolarized (P = 0), fully left-handed polarized (P = —1), and fully 
right-handed polarized (P = +1). Since we give explicit analytical expres- 
sions, suitable modification to more realistic polarizations would be straight- 
forward. 

The results for the case of polar distributions of the decay lepton, re- 
stricted to t/s values of 400 and 800 GeV, were reported earlier in a short 



paper We present here more details, as well as azimuthal distributions. 
We also consider the possibility of a higher energy option of the linear col- 
lider, with v/i = 1500 GeV. 

Our main result may be summarized as follows. By and large the dis- 
tribution in the polar angle Q\ of the secondary lepton w.r.t. the e~ beam 
direction is unchanged in shape on inclusion of QCD corrections in SGA. 
The Q\ distribution for y/s = 400 GeV is very accurately described by overall 
multiplication by a K factor (K = 1 + k > 1), except for extreme values of 
9i, and that too for the case of P = +1. For the other energies considered, « 
continues to be slowly varying function of Q\. This has the important conse- 
quence that earlier results on the sensitivity of lepton angular distributions or 
asymmetries to anomalous top couplings, obtained for y/s values around 400 
GeV without QCD corrections being taken into account, would go through 
by a simple modification by a factor of 1/ \f~K . The triple distributions 
in 9 t , Qi and (f>i show an asymmetry around 0/ = 180°, which is not present 
at Born level. 
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2 Derivation of angular distributions to 0(a s ) 



We first obtain expressions for helicity amplitudes for 

e~(p e -) + e + (p e +) -> t(p t )+t(pt) 



going through virtual 7 and Z in the e + e~ cm. frame, including QCD 
corrections in SGA. The starting point is the QCD-modified jit and Zit 
vertices obtained earlier (see, for example, pi]]). We can write them [2J5 
in the limit of vanishing electron mass as 



r 7 = e 
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The form factors A and i3 are given to order a s in SGA by 
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Hen = a, — - — log ■ 
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where a s = a s /(3ii), (3 = \Jl — Amf/s, and Li 2 is the Spence function. ReA 
in eq. (0) contains the effective form factor for a cut-off u; max on the gluon en- 
ergy after the infrared singularities have been cancelled between the virtual- 
and soft-gluon contributions in the on-shell renormalization scheme. Only 
the real part of the form factor A has been given, because the contribution 
of the imaginary part is proportional to the Z width, and hence negligibly 
small 



, |29| . The imaginary part of B, however, contributes to the triple 
angular distribution. 

The vertices in eqs. (fj) and (||) can be used to obtain helicity amplitudes 
for e + e~ — > tt, including the contribution of s-channel 7 and Z exchanges. 
The result is, in a notation where the subscripts of M denote the signs of the 
helicities of e~, e + , t and t, in that order, 
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where 8 t is the angle the top-quark momentum makes with the e~ momen- 
tum, 7 = 1/vT — /3 2 , and tl,r are related to the left- and right-handed Zee 
couplings, and are given by 
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Since we are interested in lepton distributions arising from top decay, 
we also evaluate the helicity amplitudes for t — > bl + ui ( or t — > bl~Vi), 
which will be combined with the production amplitudes in the narrow-width 
approximation for t (i) and W + (W~). In principle, QCD corrections should 
be included also in the decay process. However, in SGA, these could be 
written in terms of effective form factors [32]. As found earlier [15], |1|] , 
in the linear approximation, these form factors do not affect the charged- 
lepton angular distribution. Hence we need not calculate these form factors. 
Nevertheless, to illustrate this point we will discuss an arbitrary tbW vertex, 
in the limit of vanishing bottom mass. 

The most general tbW vertex may be written as 



1 t^bW 



^V tb u(p b ) [YUilPl + HrPr) 
" -a^( Pt - Pb Uf 2L P L + f 2R P R ) 



m w 



u(p. 



tj- 



(19) 



where Pl,r — |(1 T7s), an d V tb is the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix element, which we will take to be equal to one. jm and jit, do not 
contribute in the limit of vanishing b mass, and can be omitted. Moreover, 
we will set fn = 1, the value at tree level. The result does not depend on 
the value of fiL, as we will see soon. The decay helicity amplitudes in the t 
rest frame can be found in 



16 



and we do not repeat them here. 
It was found in Jl6| that, in this case, the spin density matrix for t — > 
bl + ui, on integration over the 6-quark azimuthal angle, is the same as that for 
the Born vertex, except for an overall factor (l + Ref 2 R ) multiplying all 
the matrix elements. Moreover, after performing a boost to the lab. frame, 
and integrating over the charged-lepton energy, that the f 2 R dependent factor 
cancelled against the corresponding factor in the correction to the top width. 
If we had not set f iL = 1, the fn dependence would also have likewise 
cancelled. 

The final result for the angular distribution in the lab. frame can be 
written as 
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- /3 2 ) sin 9 t sin 9 X sin . 



(20) 



where (9 t and (9; are polar angles of respectively of the t and Z + momenta 
with respect to the e~ beam direction chosen as the z axis, and <pi is the 
azimuthal angle of the l + momentum relative to an axis chosen in the it 
production plane. B\ is the leptonic branching ratio of the top. Q t i is the 
angle between the t and Z + directions, given by 
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Integrating over 0/ and # t ; we get the Q\ distribution: 
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Figure 1: The distribution in #/ with and without QCD corrections for (a) y/s = 
400 GeV, (b) </s = 800 GeV and (c) ^ = 1500 GeV plotted against 9 h for e 
beam polarizations P = 0, — 1, +1 in each case. 
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3 Numerical Results 



After having obtained analytic expressions for angular distributions, we now 
examine the numerical values of the QCD corrections. 

We use the parameters a = 1/128, a s {m 2 z ) = 0.118, mz = 91.187 GeV, 
m w = 80.41 GeV, m t = 175 GeV and sin 2 W = 0.2315. We consider 
leptonic decays into one specific channel (electrons or muons or tau leptons), 
corresponding to a branching ratio of 1/9. We have used, following [p9fl . 
a gluon energy cut-off of cu max = (y/s — 2m t )/5. While qualitative results 
would be insensitive, exact quantitative results would of course depend on 
the choice of cut-off. 

In Figs. 1 (a)-(c) we show the single differential cross section d c ^ in 
picobarns with and without QCD corrections, for three values of ^/s, viz., 
(a) 400 GeV, (b) 800 GeV and (c) 1500 GeV, and for e" beam polarizations 
P = 0,-1, +1. It can be seen that the distribution with QCD corrections 
follows, in general, the shape of the lowest order distribution. 

In Figs. 2 (a)-(c) is displayed the fractional deviation of the QCD- 
corrected distribution from the lowest order distribution: 

(Q\ _ ( d&Born\ ( dlTsGA _ do Born\ 

M l) ~ [dcosej \dcos6t ~ dcosej ' 1 j 

It can be seen that k(9i) is independent of 9i to a fair degree of accuracy for 
v/i = 400 GeV. 

In Figs. 3 (a)-(c) we show the fractional QCD contributions (-Fsga — 
-^Bom) /inborn where F(9i), is the normalized distribution: 

FW = I(^_). (28) 
a \d cos 6 1 J 

It can be seen that the fractional change in the normalized distributions for 
V^i = 400 GeV is at most of the order of 1 or 2% (except in the case of 
P = +1, for $i > 160°). For the other values of ^/s, it is even smaller. This 
implies that QCD corrected angular distribution is well approximated, at the 
per cent level, by a constant rescaling by a K factor. 

We now discuss the triple differential cross section, which corresponds 
to a distribution in 9t, Q\ and (pi, for which we have obtained an analytic 



expression, eq. (p0|) . This distribution, unlike the 8[ distribution, requires 
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Figure 2: The fractional QCD contribution k{6i) denned in the text for (a) y/s = 
400 GeV, (b) y/s = 800 GeV and (c) yfi = 1500 GeV plotted as a function of h 
for P = 0,-1, +1. 
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Figure 3: The fractional QCD contribution in normalized angular distributions, 
F(9i) defined in the text, for (a) y/s = 400 GeV, (b) y/s = 800 GeV and (c) 
y/s = 1500 GeV plotted as a function of 9i, for P = 0, —1, +1. 
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the knowledge of the direction of t. In fact, without knowing the direction 
of the top momentum, it is not even possible to define (pi. However, the 
magnitude of the momentum need not be determined precisely. 

We have plotted in Fig. 4(a) d 3 a j(d cos 9 t d cos 9id(pi), including QCD 
correction, against (pi for different values of 9i, for a fixed value of 9 t = 30°, 
for t/s = 400 GeV, and for P = 0. Fig. 4(b) gives the corresponding 
fractional deviations from the lowest-order values, 

d 3 a SGA d 3 a B 

orn 

d cos 9 t d cos 9 \d(p x d cos 9 t d cos 9id(pi 

(29) 

It can be seen that the QCD contribution remains around 28% for various 
values of Q\ and (pi. The triple distributions show an asymmetry around 
(pi = 180°, which is not present at Born level. This asymmetry is not visible 
in the Fig. 4 (a), but is very clear in Fig. 4 (b), where fractional QCD 
contributions to the normalized distributions are plotted. 

The same plots as in Fig. 4 (a), (b) are repeated for y/s = 800 GeV in 
Figs. 5 (a), (b) and for = 1500 GeV in Figs. 6 (a), (b). 

We have not shown the plots for other values of 9 t because they do not 
illustrate anything new. We have also not discussed triple distributions for 
the case of polarized beams for the sake of brevity. 

4 Conclusions 

To conclude, we have obtained in this paper analytic expressions for angular 
distributions of leptons from top decay in e + e~ — > ti, in the e + e~ cm. frame, 
including QCD corrections to order a s in the soft-gluon approximation. The 
distributions are in a form which can be compared directly with experiment. 
In particular, the single differential 6i distribution needs neither the recon- 
struction of the top momentum direction nor the top rest frame. The triple 
differential distribution does need the top direction to be reconstructed for 
the definition of the angles. However, the magnitude of the top momentum 
need not be determined. In either case the results do not depend on the 
choice of spin quantization axis. 

We find that the 9i distributions are well described by rescaling the zeroth 
order distributions by a factor K which for ^Js = 400 GeV is roughly inde- 
pendent of 9i, except for extreme values of for the case of right-handed 
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Figure 4: (a) The triple differential distribution including the QCD contri- 
bution for y/s = 400 GeV and 9 t = 30°. (b) The fractional QCD contribu- 
tion K,(9 t , 9i, 4>i) to the triple differential distribution for y/s = 400 GeV and 
9 t = 30°. 
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Figure 5: (a) The triple differential distribution including the QCD contri- 
bution for y/s = 800 GeV and 9 t = 30°. (b) The fractional QCD contribu- 
tion n(9t, 0i, 4>i) to the triple differential distribution for y/s = 800 GeV and 
9 t = 30°. 
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9 t = 30°. 
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polarized electron beam. For other values of y/s, it is a slowly varying func- 
tion of Q[. 

The triple distributions in 6 t , 9i and (pi show an asymmetry about (pi = 
180°, which is absent at tree level, and would be a good test of the QCD 
corrections. 

It would be useful to carry out the hard-gluon corrections explicitly and 
check if the soft-gluon approximation used here gives correct quantitative 
results. 
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